Search Results

The objective of this paper is to assess the effect of snow density on tire-snow interaction in the presence of uncertainty. The snow-depth dependent finite element analysis (FEA) and semi-analytical models we have developed recently can predict tire-snow interfacial forces at a given density under combined slip conditions. One drawback of the models is that they are only applicable for fresh, low-density snow due to the unavailability of a density-dependent snow model. In reality, the snow density on the ground can vary between that of fresh snow to heavily compacted snow that is similar to ice. Even for fresh snow on the ground, as a vehicle moves forward, the rear wheels experience higher snow densities than the front wheels. In addition, being a natural material, snow's physical properties vary significantly even for the same density.

Existing methods for design optimization under uncertainty assume that a high level of information is available, typically in the form of data. In reality, however, insufficient data prevents correct inference of probability distributions, membership functions, or interval ranges. In this article we use an engine design example to show that optimal design decisions and reliability estimations depend strongly on uncertainty characterization. We contrast the reliability-based optimal designs to the ones obtained using worst-case optimization, and ask the question of how to obtain non-conservative designs with incomplete uncertainty information. We propose an answer to this question through the use of Bayesian statistics. We estimate the truck's engine reliability based only on available samples, and demonstrate that the accuracy of our estimates increases as more samples become available.

An accurate and efficient Monte Carlo simulation (MCS) method is developed in this paper for limit state-based reliability analysis, especially at system levels, by using a response surface approximation of the failure indicator function. The Moving Least Squares (MLS) method is used to construct the response surface of the indicator function, along with an Optimum Symmetric Latin Hypercube (OSLH) as the sampling technique. Similar to MCS, the proposed method can easily handle implicit, highly nonlinear limit-state functions, with variables of any statistical distributions and correlations. However, the efficiency of MCS can be greatly improved. The method appears to be particularly efficient for multiple limit state and multiple design point problem. A mathematical example and a practical example are used to highlight the superior accuracy and efficiency of the proposed method over traditional reliability methods.

An efficient simulation-based method is proposed for the reliability analysis of a vehicle body-door subsystem with respect to an important quality issue -- wind noise. A nonlinear seal model is constructed for the automotive wind noise problem and the limit state function is evaluated using finite element analysis. Existing analytical as well as simulation-based methods are used to solve this problem. A multi-modal adaptive importance sampling method is then developed for reliability analysis at system level. It is demonstrated through this industrial application problem that the multi-modal adaptive importance sampling method is superior to existing methods in terms of efficiency and accuracy. The method can easily handle implicit limit-state functions, with variables of any statistical distributions.

This paper presents the development of surrogate models (metamodels) for evaluating the bearing performance in an internal combustion engine. The metamodels are employed for performing probabilistic analyses for the engine bearings. The metamodels are developed based on results from a simulation solver computed at a limited number of sample points, which sample the design space. An integrated system-level engine simulation model, consisting of a flexible crankshaft dynamics model and a flexible engine block model connected by a detailed hydrodynamic lubrication model, is employed in this paper for generating information necessary to construct the metamodels. An optimal symmetric latin hypercube algorithm is utilized for identifying the sampling points based on the number and the range of the variables that are considered to vary in the design space.

A fast and accurate journal bearing elastohydrodynamic analysis is presented based on a finite difference formulation. The governing equations for the oil film pressure, stiffness and damping are solved using a finite difference approach. The oil film domain is discretized using a rectangular two-dimensional finite difference mesh. In this new formulation, it is not necessary to generate a global fluidity matrix similar to a finite element based solution. The finite difference equations are solved using a successive over relaxation (SOR) algorithm. The concept of “Influence Zone,” for computing the dynamic characteristics is introduced. The SOR algorithm and the “Influence Zone” concept significantly improve the computational efficiency without loss of accuracy. The new algorithms are validated with numerical results from the literature and their numerical efficiency is demonstrated.

A comprehensive formulation is presented for the dynamics of a rotating flexible crankshaft coupled with the dynamics of an engine block through a finite difference elastohydrodynamic main bearing lubrication algorithm. The coupling is based on detailed equilibrium conditions at the bearings. The component mode synthesis is employed for modeling the crankshaft and block dynamic behavior. A specialized algorithm for coupling the rigid and flexible body dynamics of the crankshaft within the framework of the component mode synthesis has been developed. A finite difference lubrication algorithm is used for computing the oil film elastohydrodynamic characteristics. A computationally accurate and efficient mapping algorithm has been developed for transferring information between a high - density computational grid for the elastohydrodynamic bearing solver and a low - density structural grid utilized in computing the crankshaft and block structural dynamic response.

High vibration levels at the rear bearing cap and oil pump were observed in dyno tests for a particular design of a V6 engine at a rated speed of 4800 r/min. It was found experimentally that the crankshaft-flywheel assembly had a bending resonance at 240 Hz which was excited at around 4800 r/min by 3rd order forces on the crankshaft. A newly developed crankshaft system model (CRANKSYM) was used to analytically verify the above finding and propose possible solutions to the problem. CRANKSYM can perform a coupled analysis among the crankshaft structural dynamics, main bearing hydrodynamics and engine block flexibility. It considers the flywheel dynamics (including the gyroscopic effect), belt loads, crankshaft “bent” and block misboring, and the anisotropy of the block flexibility as seen from a rotating crankshaft. It can also calculate the dynamic stresses on the crankshaft throughout the whole engine cycle. A brief description of CRANKSYM is given in the paper.

The side-load capacity of an unlubricated ringless piston assembly was calculated for a four-stroke I.C. engine for possible application to a Low-Heat-Rejection engine. A compressible flow analysis for an ideal gas was used to calculate the flow at the piston-cylinder interface. The analysis accounts for the gas inertia effect and possible choked flow. The piston side-load capacity is compared with the actual side thrust on the piston, generated through the combination of piston inertia, cylinder pressure, and connecting rod angularity. It was found that a carefully designed ringless piston can support the side load of a slider-crank mechanism during the power stroke of a four-stroke engine. However, during the remaining strokes, the load can be supported only partially.

A theoretical analysis is presented for calculating the trajectory of a ringless piston within the cylinder clearance of an I.C. engine with a crosshead design. The flexible polytope unconstrained minimization method is used to find the equilibrium between the piston rod structural elasticity and the gas-film hydrodynamics at the piston-cylinder interface. It was found that even when the crosshead bearing is not concentric with the cylinder (i.e., there is an initial eccentricity between the piston and cylinder centerlines), the piston does not touch the cylinder wall during an engine cycle. However, this happens only when a carefully designed piston skirt profile and piston rod length and diameter are used.